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Abstract 

We investigate the solar neutrino problem assuming simultaneous presence of 
MSW transitions in the sun and neutrino decay on the way from sun to earth. We 
do a global x 2 -analysis of the data on total rates in CI, Ga and Superkamiokande 
(SK) experiments and the SK day-night spectrum data and determine the changes 
in the allowed region in the Am 2 — tan 2 9 plane in presence of decay. We also 
discuss the implications for unstable neutrinos in the SNO experiment. 
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1 Introduction 



The global analysis of the total rates measured in the CI, Ga |lj and SK experiments 
and the day-night spectrum data of SK || indicate that the large mixing angle 
MSW solution gives the best description of the solar neutrino data 0. However 
before a particular solution can be established one should rule out other possibili- 
ties. In this spirit people have considered various non-standard neutrino properties 
and their implication for the solar neutrino problem [[| . In this paper we consider 
a scenario where neutrinos are allowed to decay on their way from sun to earth 
after undergoing MSW transitions in the sun. Such a possibility was discussed 
earlier in |J. 

We consider two flavor mixing between v e and v^j v T with the mass eigenstates 
v\ and U2- We assume that the heavier mass state V2 is unstable, while the lighter 
neutrino mass state v\ has lifetime much greater than the sun-earth transit time 
and hence can be taken as stable. There are two possible non-radiative decay 
modesQ 

• Model 1: If neutrinos are Dirac particles one has the decay channel v-i 
vm+(f), where v\r is a right handed singlet and <fi is an iso-singlet scaler. Thus 
all the final state particles for this model are sterile and there is no distinct 
signature of this decay apart from in disappearance experiments. This model 
is discussed in |7]]. In this model a light scalar boson with lepton number 
-2 and a singlet right handed neutrino is added to the standard model. The 
neutrino coupling to this scalar boson is given by g2i^R 1 C~ 1 i'R 2 , C being the 
charge conjugation operator. 

• Model 2: If neutrinos are Majorana particles, the decay mode is z/ 2 — > V\ + J, 
where J is a Majoron, produced as a result of spontaneous breaking of a global 
[/(l)i e _L symmetry M. In this model the neutrino masses are generated by 
extending the higgs sector of the standard model. Though the original triplet 
majoron model proposed by Gelmini and Roncadeli H is ruled out from the 



LEP data on Z decay to invisible modes p0| , the model discussed in jS|, 
which needs two additional triplet and one singlet scalar boson in the theory, 
can avoid the conflict with the LEP data, and at the same time predict fast 
enough neutrino decay necessary to solve the solar neutrino problem. In this 
model the v\ can be observed with a probability \U e i\ 2 and as a v^jvr 

with a probability |{7 e2 | 
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1 Radiative decays are severely constrained! 
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In both the decay scenarios the rest frame lifetime of v-i is given by [|TT] 



^_ _ 167rm 2 (l + m 1 /m 2 ) 2 ^ 
g 2 Am 2 

where g is the coupling constant, m, is the U{ mass and Am 2 = m 2 ,— m\ . Assuming 
vri2 » mi the equation ([!]) can be written as 

g 2 Am 2 ~ lQira (2) 

where a is the decay constant related to tq as a = 7712/70. If we now incorporate 
the bound g 2 < 4.5 x 10~ 5 as obtained from K decay modes [O] we get the bound 



Am 2 > 10 6 a 0. Since for a typical neutrino energy of 10 MeV, one starts getting 
decay effects over the sun-earth distance for a ~ 10~ 13 eV 2 , the corresponding limit 
on Am 2 is Am 2 ~ 10~ 7 eV 2 . In an earlier paper we had considered high values of 
Am 2 (> 10~ 3 eV 2 ) so that the matter effects inside the sun can be neglected and 
one can have decay as well as Am 2 independent average oscillations [14|. In this 
paper we consider 10 -6 < Am 2 < 10~ 3 eV 2 such that the matter effects inside the 
sun are important. We incorporate the earth matter effects as well. 

In section 2 we discuss the usual two flavor MSW solutions to the solar neutrino 
problem. We perform a x 2 -analysis to the global data on rates and SK day-night 
spectrum and present the allowed regions in the Am 2 — tan 2 9 plane. In section 3 
we introduce the possibility of having one of the neutrino states, v% to be unstable. 
We look for the effects of decay on the MSW solutions and present the allowed areas 
at various values of the decay constant. We show that for the majoron decay model 
(model 2) the SK data on v e — p events restricts the allowed values of the decay 
constant to extremely low values. In section 4 we discuss the implications of non- 
zero neutrino decay for the SNO experiment. We finally present our conclusions 
in section 5. 



2 MSW effect for stable neutrinos 

As was first indicated in JT3 , interaction of electron neutrinos propagating through 



the sun modifies the vacuum mixing angle 9 to matter mixing angle 9m where, 

Am 2 sin 29 , x 

tan2# M = = . (3) 

Am 2 cos 29 -2V2G F n e E v ; 



5 Since the m ee element of the mass matrix is zero in the majoron decay model considered, 
the 0^-majoron (3j3 decay does not take place in this model Q and the more stringent bound on 
g from majoron emission in (3(3 decay p3] is not applicable. 
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Here n e is the ambient electron density, E the neutrino energy, and Am 2 (= 
wij — m\) the mass squared difference in vacuum. The vanishing of the denomi- 
nator in eq. ([5]) defines the resonance condition, where the matter mixing angle is 
maximal. 

The probability amplitude of survival for an electron neutrino is given by, 

A - A® A vac A® -l- A Q A vac A® (A\ 

where Af k gives the probability amplitude of u e —>■ transition at the solar surface, 



A0 



A% = a%e~^ (5) 

a® 2 = Pj sin 2 9 M + (1 - Pj) cos 2 6 M = 1 - a© 2 (6) 

Pj is the non-adiabatic level jumping probability between the two mass eigenstates 
for which we use the standard expression from fT6| . 

exp(-7sin 2 ^) -exp(-7) 

Pj = ; t — x { i ) 

1 - exp(-7) 



7 = 7T- 



Am 2 



d lnn e 
dr 



A11 C is the transition amplitude from the solar surface to the earth surface, 

a vac _ -iE k (L-R ) 
A kk — e 



(9) 



where L is the sun-earth distance and R & the radius of the sun. A® e denotes the 
Vk —> v e transition amplitudes inside the earth. We evaluate these amplitudes 
by assuming the earth to consist of two constant density slabs. The v e survival 
probability is given by 



a 



©2, 



el 



A® 



+ a e2 l /1 2el 



+2af ia %Re[Af e Al*e 



i{E i -E 1 )(L-R @ ) e i(<h,<3-<t>i,e)] 



(10) 



It can be shown that the square bracketed term containing the phases average out 
to zero in the range of Am 2 in which we are interested ||17|| . 

Using eq. fllPf) as the probability we have done a x 2 -analysis of the current solar 
neutrino data on total rates from CI, Ga[] and SK and the 18+18 bins of data on 

6 We use the weighted average of the rates from Sage, Gallex and GNO. 
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the day-night electron energy spectrum from SK||18||. The details of the code used 



can be found in [19|. We use the 1117 days data of SK and incorporate the BBPOO 
solar model The theory errors and their correlations in the analysis of total 

El ' 



rates is included as in 



In addition to the astrophysical uncertainties included 
in |]19| for the analysis of the total rates, we have included the uncertainty in the 
So-factor for the reaction w O(p, , ~f) 17 F | 2~0[| . For the day-night spectrum analysis 



we have included the correlation between the systematic errors of the day and 
the night bins. As in [22| we vary the normalisation of the spectrum as a free 
parameter which avoids the overcounting of the rates and spectrum data for SK. 
Hence for the day-night spectrum analysis we have (36 - 1) degrees of freedom 
(d.o.f) while for the total rates we have 3, which makes a total of 38 d.o.f for the 
rates+spectrum analysis with no oscillation. For MSW analysis with two active 
neutrino flavors we present in Table 1 two sets of results 

• with the 8 B normalisation factor Xb in the total rates held at the SSM value. 

• with the 8 B normalisation factor Xb in the total rates kept as a free param- 
eter. 





Nature of 
Solution 


Am 2 
in eV 2 


tan 2 6» 


V 2 . 

Amin 


Goodness 
of fit 


Xb fixed 
at SSM value 


SMA 


5.48 x KT 6 


5.79 x 10" 4 


43.22 


19.01% 


LMA 


4.17 x 1(T 5 


0.35 


37.33 


40.78% 


LOW 


1.51 x l(T y 


0.64 


39.54 


31.48% 


X B varying 


SMA 


5.35 x 10" 6 


4.34 x 10" 4 


37.98 


33.51% 


LMA 


4.22 x 1(T 5 


0.25 


34.20 


50.67% 


LOW 


1.51 x l(T y 


0.64 


39.88 


26.20% 



Table 1: The best-fit values of the parameters, Xmim an d the goodness of fit 
from the global analysis of rates and the day-night spectrum data for MSW 
analysis involving two neutrino flavors. 



Thus the best-fit comes in the LMA region. In fig. 1 we plot the 90%, 95% and 
99% CL. allowed regions for the two-flavor MSW transition to an active neutrino. 
All the contours have been drawn with respect to the global minimum with the 
Xb fixed at the SSM value. We have also done a x 2 " anar ysis f° r v e — v s MSW 
conversion, v s being a sterile neutrino, for both Xb fixed at SSM and Xb varying 
freely. The best-fit values of parameters, Xmin anc ^ ^ e goodness of fit (g.o.f) are 

• Am 2 = 3.74 x 10~ 6 eV 2 , tan 2 # = 5.2 x 10" 4 , xLn = 44 - 85 > g-°- f = 
14.79%, X B fixed at SSM value. 



Am 2 = 3.71 x 10~ 6 eV 2 , tan 2 9 = 4.72 x 10 
15.53%, X B free. 



v 2 . 

A,min 



= 43.42, g.o.f = 
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3 MSW effect and unstable neutrinos 



If a neutrino of energy E decays while traversing a distance L then the decay term 
exp(— aL/E) gives the fraction of neutrinos that survive, where a is the decay 
constant discussed before. In fig. 2 we plot this exp(— aL/E) as a function of 
a for two different energy values. For a very small the exponential term is ~ 1 
and there is no decay. As a increases one starts getting decay over the sun-earth 
distance only for a ~ 10~ 13 eV 2 . As we have discussed in the introduction, since 
from K-decay Am 2 > 10 6 a, one can have simultaneous MSW and decay for 10~ 7 
eV 2 < Am 2 < 10~ 3 eV 2 . Below this Am 2 the bound on the coupling constant 
coming from K-decay restricts the decay constant a to be small enough so that 
decay effects are negligible over the sun-earth distance. So henceforth we will be 
concerned with only the LMA and SMA solutions, the LOW region remaining 
unaffected due to decay. For Am 2 > eV 2 there will not be any MSW effect. 

The probability amplitude for v e survival in presence of neutrino decay is again 
given by eq.(f|), with Af k , A® e and A™ c as before, the only change being for A^ c 
since the v 2 decays on its way from the sun to earth. is given by 

A vac = e - iE 2(L-R@) e -a(L-R @ )/E 2 

Then, the v e survival probability is given by (ignoring the phase part), 



P ee = a^\Aff + a% 2 \AZ\ 2 e-^ L - R ^ (12) 

The day-time probability (i.e. without the earth effect) is given by 

P e d e ay = cos 2 9[Pj sin 2 9 M + (1 - Pj) cos 2 9 M ] 

+ sin 2 0[(1 - Pj) sin 2 9 M + cos 2 9 M Pj]e~ 2a(L - R ^ (13) 

From eq. (O) we note that the the decay term appears with a sin 2 9 and is therefore 
appreciable only for large enough 9. Thus we expect the effect of decay to be 
maximum in the LMA region. This can also be understood as follows. The v e are 
produced mostly as z/ 2 in the solar core. In the LMA region the neutrinos move 
adiabatically through the sun and emerge as v 2 which eventually decays. For the 
SMA region on the other hand Pj is non-zero and u e produced as v 2 cross over 
to V\ at the resonance and come out as a V\ from the solar surface. Since V\ is 
stable, decay does not affect this region. Including the earth effect the survival 
probability can be expressed as 

day (sin 2 9 - P 2e ) (P^ + e -ML-R Q )/E 2 (p da y _ 1)} 

ee cos 2Q _ stn 2 6e ~2 a (L-R Q )/E 2 ^ > 
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P2e is the probability of the second mass eigenstate to get converted into the v e 
state at the detector. 

In fig. 3 the solid lines show the survival probability versus neutrino energy 
for three different values of a with Am 2 and tan 2 9 fixed at the best fit value of 
the LMA solution. The introduction of decay reduces the survival probability as 
the second term in eq.(|T2D reduces with a. Decay also brings about a deviation 
from a flat probability towards the high energy end. The dotted line gives the 
probability for the same Am 2 but a higher value of tan 2 9 with a = 10 -11 eV 2 . We 
observe that there is a huge drop in the survival probability initially, followed by a 
sharp increase with energy. A comparison between the different curves shows that 
inclusion of decay results in the reduction of the survival probability as well as an 
energy distortion and both these effects increase with the mixing angle 9. 

We next perform a x 2 -analysis of the total rates and the day-night spectrum 
data for non-zero decay. The best-fit comes for a = 0, corresponding to no decay 
with the Am 2 and tan 2 9 as given in Table 1. However non-zero values of a giving 
finite decay probability also gives acceptable fit. For small values of 9 as in the 
SMA region, the fraction of z/2 in the solar neutrino beam is very small. As a result 
very few neutrinos decay and all values of a may be allowed. On the other hand in 
the LMA region the 9 is high and large number of neutrinos can decay producing a 
distortion in the energy spectrum. Therefore the data can put some restrictions on 
the allowed values of a. In fig. 4 we plot A% 2 (= x 2 — Xmin) vs a keeping the other 
two parameters free (in the LMA regime). We also show the 99% C.L. limit for 
three parameters by the dotted line. From fig. 4 we see that in the LMA region, 
values of a upto 3.5 x 1CT 11 eV 2 are allowed at 99% C.L. from the global analysis. 

In fig. 5 we show the allowed regions in the Am 2 — tan 2 9 plane for various 
fixed values of a. The first panel is for a = which is the case of no decay and 
the contours are the same as those presented in fig. 1, modulo the difference in the 
definition of the C.L. as now there are three parameters whereas in fig. 1 we had 
two. The other three panels are for different non-zero allowed values of a obtained 
from fig. 4. 

As expected, we find that the SMA allowed regions remain mostly unaffected 
due to a non-zero a. The LMA solution on the other hand is appreciably affected 
and the allowed area shrinks as a increases. This effect is seen to be more pro- 
nounced for larger tan 2 9. The reason for this can be traced down to the fact that 
decay results in distortion of the high energy end of the neutrino spectrum and 
this distortion is more for higher tan 2 9, a feature explicit in fig. 3. As a result 
higher values of tan 2 9 are disfavored by the global data as a increases. 

It is to be noted that although in the LOW region the sin 2 9 is high, the Am 2 is 
small and appreciable decay over the sun-earth distance is not obtained if one has 
to be consistent with the bounds on the coupling constant from K-decays. Because 
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of this the LOW region is not plotted in fig. 5. 

The results presented in this section are in general valid for both the decay 
models though in model 2 the v 2 decays to a i>\, which can interact with the 
electrons in SK detector as u e and v x . But since the V\ is degraded in energy 0. 
the effect of these additional {y e — e) and {p x — e) scattering events in SK is not 



significant [|14| and the final results remain the same. However additional constraint 
on model 2 come from the v e p — > e + n events, which contribute to the background 
in SK. From the absence of a significant contribution above the background, the 
bound obtained on the total flux of i> e from 8 B neutrinos is <lv e ( 8 -B) < 1.8 x 10 5 
cm~ 2 s _1 which translates to a bound on the probability P eS < 3.5% ||23|| . In fig. 6 



we plot the conversion probability P eS vs a. For each a we find the Xmin an d plot 
P eS for the corresponding the best fit Am 2 and tan 2 9. We also show the allowed 
limit of Pee = 0.035. So from this constraint for model 2 only a < 5.8 x 10~ 13 eV 2 
remain allowed. 



4 Implications for SNO 

The main detection processes in the heavy water of SNO are 

h> e + d— > p + p + e~ (15) 



v x + d^p + n + v x (16) 

The first is a charged current (CC) reaction with energy threshold of 1.44 MeV 
and the second one is a neutral current (NC) process with an energy threshold 
of 2.23 MeV. While the CC is sensitive to only u e , the NC process is sensitive to 
neutrinos and antineutrinos of all flavors. The rate of iv e ~ d) CC events recorded 
in the detector is given by 

_ fdE u \ Ue (E„)a C c(E u )(Pee) 
CC J dE v KXEu)(Tcc{E v ) 1 ' 



acc = [ dE A r dE T R(E A , E T ) f dE v ^§^> (18) 

JE Ath JO J CLEt 

where X Ue is the normalised 8 B neutrino spectrum, (P ee ) is the time averaged 
v e survival probability, da^d/dEr is the differential cross section of the [y e — d) 
interaction f24|, Et is the true and E A the apparent (measured) kinetic energy of 
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the recoil electrons, E& th is the detector threshold energy which we take as 5 MeV 
and K(Ea,Et) is the energy resolution function which is assumed to be Gaussian 



^(0.348^) 



1 



exp 



( 



0.2A2E T MeV 



(Et — E A ) 2 



) 



(19) 



The NC event rate is given by 



Rnc = 



J dE u \ Ue (E u )a NC (E u )V 
J dE v \ Ve {E v )a NC {E v ) 



(20) 



where V = (P ee ) + (Pep) f° r two flavor v e — oscillation. The left hand panels 
in fig. 7 give the range of predicted values of Rcc, Rnc and the double ratio 
Rnc/Rcc at 99% C.L., from the global MSW analysis of the solar neutrino data. 
The black dots in these figures represent the expected rates in SNO for the local 
best fit values of the parameters obtained. As the neutral current interaction is 
flavor blind, the Rnc remains 1 for oscillations to active neutrinos for which V — 1. 
But for oscillations to sterile neutrinos V = (P ee ) and Rnc = 0.78 for the best fit 
case. 

The right hand panels in fig. 7 show the the corresponding rates in SNO for the 
LMA region if the neutrinos are assumed to be unstable (the SMA solution remains 
largely unaffected as discussed in the previous section). We consider the decay 
model 1 and present the 99% C.L. predicted range of values of the rates for three 
different allowed values of the parameter a, along with the case for a — 0, which 
is the same as the LMA region of global MSW analysis, modulo the difference 
in the definition of the C.L.. Also shown are the rates for the best fit value of 
Am 2 and tan 2 9 for various fixed values of a. For each fixed value of Am 2 and 
tan 2 # one expects Rcc to fall sharply with alpha. But since the best fit of the 
global analysis for increasing a corresponds in general to a higher Am 2 , that raises 
the value of Rcc- Hence Rcc m SNO alone may not be able to distinguish the 
decay scenario from MSW conversion to active stable neutrinos. The value of Rnc 
suffers a marked decrease from 1 as the decay constant a increases, since for the 
decay model 1 the v e decays to sterile particles which cannot be detected in SNO. 

From fig. 7 we see that if Rnc is below 1, one can have either MSW conversion 
to sterile neutrinos, or the possibility of unstable neutrinos. But the ambiguity 
between the last two cases can be lifted if one next looks at the value of Rcc, 
which is much lower for the case of unstable neutrinos than what one would get 
for MSW conversion to sterile neutrinos. We also display the expected range of 
Rnc/Rcc f°r various fixed values of the decay constant. It is clear from the figure 
that from the value of Rnc/Rcc m SNO one can clearly distinguish the case of 
decay from the case of MSW transition to sterile neutrinos. 
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In the decay model 2 the v e decay to i/j, which can in principle show up both 
in the NC events in SNO as well as in the charged current absorption reaction 

v e + d->n + n + e + (21) 

But since the absence of (u e — p) events in SK restrict a < 5.8 x 10~ 13 eV 2 , both 
the additional contribution to the NC event as well as the {p e — d) absorption is 
negligible in SNO. 



5 Conclusions 

We have investigated the effect of neutrino decay on the MSW solution to the solar 
neutrino problem. There are two factors which control the effect of decay 

• The mixing angle 9 which determines the fraction of the unstable component 
in the v e beam 

• The decay constant a which determines the decay rate 

In order for the neutrinos to decay over the earth sun distance a should be > 10~ 13 
eV 2 . Then eq.(H) an d the bounds on the coupling constant from K decays restricts 
Am 2 ~ 10~ 7 eV 2 and decay does not take place in the LOW region. We therefore 
probe the SMA and the LMA regions. We find that even the SMA region is 
not affected much because the mixing angle being small very few neutrinos would 
decay. But the effect of decay on the LMA region is significant. This is borne out 
by the x 2 - ana lysi s of the global solar data on rates and SK day-night spectrum. 
We point out that although the data prefers no decay, it still cannot rule out the 
possibility of unstable neutrinos completely. We put limits on the allowed range of 
the decay constant a and present the allowed areas in the Am 2 — tan 2 9 parameter 
space for various allowed values of a. The LMA allowed zone is seen to be severely 
constrained by decay. 

From fig. 4, values of a < 3.5 x 10~ n eV 2 are allowed, implying a neutrino rest 
frame lifetime Tq > 2 x 10~ 5 sec. Thus one may encounter decay before neutrino 
decoupling in the early universe. This may result in increasing the number of light 
neutrinos N v to greater than 3. However, the upper limit on N v can be as large as 



6 [25|. Hence our decay model is consistent with the bounds from early universe. 

From the fact that the neutrinos from SN1987A have not decayed on their way 
one gets a lower bound on the electron neutrino lifetime as r > 5.7 x 10 5 (m jye /eV) 
sec. However if one includes neutrino mixing then shorter lifetimes are allowed 



provided \U e 2\ < 0.9 fl26| . This again is consistent with our analysis. 

In this paper we assume one of the neutrino states to be unstable in vacuum 
and the values of the decay constant a are such that one has decay over earth-sun 
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distance. For these values of a decay inside the sun or the earth is negligible. 
However as was shown in [27] matter can induce neutrino decay with majoron 
emission even when neutrinos are stable in vacuum. It was shown later in [^] that 
in the presence of only standard interactions, the matter induced decay cannot 
provide a solution to the solar neutrino problem due to strong bounds on neutrino- 
majoron couplings. As we have discussed in this paper, decay models where the v 2 
decays to a v\ and a majoron is independently disfavored from the non-observation 
of (i> e — p) events in SK. 

We have looked at the implications of unstable neutrinos for the SNO detector. 
We first present the values of the CC and NC rates which SNO is expected to 
observe for MSW transitions with stable neutrinos. We have given the values for 
transitions to both active as well as sterile species. For unstable neutrinos the 
NC rate is less than 1 and comparable to Rnc f° r standard MSW SMA transition 
to sterile neutrinos. But even though the value of Rnc ma y be nearly same for 
the two scenarios, comparing the values of Rcc or Rnc/Rcc h ma Y be possible 
in principle to distinguish the decay scenario from the MSW transition to sterile 
neutrinos. 

The authors wish to thank Carlos Pena-Garay for many useful discussions during 
the development of the solar code, Y. Suzuki for sending them the SK data and H.M. 
Antia for many helpful discussions. 
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Fig. 1: The 90, 95 and 99% C.L. allowed area from the global analysis 
of the total rates from CI, Ga and SK detectors and the 1117 days SK 
recoil electron spectrum at day and night, assuming MSW conversions 
to stable sequential neutrinos. 
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Fig. 2: The decay term exp(aL/E) vs a for two different neutrino 
energies. L is taken as the earth-sun distance. 
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Fig. 3: The survival probability as a function of v e energy for different 
values of a. The solid lines are for tan 2 9 = 0.35 while the dotted line 
is for tan 2 9 = 0.85. For all the curves Am 2 = 4.17 x 10~ 5 eV 2 . 
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Fig. 4: The A% 2 versus decay constant a for the global analysis of 
total rates and the SK day and night spectrum data. Also shown by 
the dotted line is the 99% CL. limit for three parameters. 
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Fig. 5: The 90, 95 and 99% CL. allowed area from the global analysis 
of the total rates and the 1117 days SK day-night spectrum data, for MSW 
conversion of unstable neutrinos at various values of the decay constant a. 
The different fixed values of a in eV 2 are indicated. 
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Fig. 6: The conversion probability of v e to P e , P eS versus decay con- 
stant a. Also shown by the dotted line is the 99% C.L. value of P eS 
allowed from the non-observance of (P e — p) events in the SK detector. 
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Fig. 7: The predicted ranges of -Rccs -Ratc an d Rnc/Rcc 111 SNO for 
the standard MSW conversion to stable neutrinos (left hand panels) and 
for unstable neutrinos (right hand panels). The dots give the values of 
the rates for the local best fit. 
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